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Abstract 

The problem of maintaining gauge invariance when truncating the two particle irreducible (2PI) 
effective action has been studied recently by several authors. Here we give a simple and very general 
derivation of the gauge dependence identities for the off-shell 2PI effective action. We consider the 
case where the gauge is fixed by an arbitrary function of the quantum gauge field, subject only to the 
restriction that the Faddeev-Popov matrix is invertable. We also study the background field gauge. We 
address the role that these identies play in solving gauge invariance problems associated with physical 
quantities calculated using a truncated on-shell 2P1 effective action. 



1 Introduction 

Considerable progress has been achieved in the study of the dynamics of quantum fields in and out of 
equilibrium. Collective effects and long range interactions are intrinsic to deconfined QCD matter in heavy 
ion collisions. At equilibrium, the hard thermal loop effective theory is the appropriate gauge invariant 
theory to describe most collective and long range effects. On the other hand, kinetic theories provide the 
best method to describe a near equilibrium situation in the small coupling constant regime. Far from 
equilibrium and/or at large coupling constant other techniques need to be developed. A promising candidate 
is the 2P1 effective action method. Unfortunately, practical calculations necessitate the use of approximate 
versions of the exact 2P1 effective action. The approximated, or truncated, 2PI effective action is simply 
a Schwinger-Dyson resummation of the two point function. Without further rcsummation for the vertex 
functions, this resummation necessarily leads to gauge dependent results for most physical quantities. 

The study of the gauge dependence of functional methods in quantum field theory has a long history. 
One of the first contributions was made by Nielsen pp who showed that the explicit dependence of the 
1PI effective potential on the gauge parameter is compensated for by the gauge parameter dependence of 
the expectation value, and that the effective potential is gauge parameter independent. In (213) a general 
functional formalism was used to derived the gauge dependence of the associated n-point functions. These 
results were then used to address the gauge fixing dependence of the one loop QCD plasma damping rate. 

In this paper we derive an expression for the gauge fixing dependence of the exact, off-shell 2PI effective 
action for any gauge fixing function that has an invertible Faddeev-Popov matrix. Using this expression for 
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the exact effective action we wiil analyze the gauge dependence of the truncated effective action and show 
that the gauge dependence always occurs at higher order, within any self-consistent truncation scheme. This 
verifies the expectation that, in complete analogy with the 1PI formalism, gauge invariance problems that 
occur in specific calculations using the 2PI formalism should ultimately be traceable to an inconsistency in 
the approximation scheme. In addition, we will show that the method we use is easily generalized to the case 
of background field gauges. The background field gauge is of interest because of the recent results of Mottola 
0] . Mottola suggested a modified form of the 2PI effective action in which Ward identities for background 
and quantum gauge fields are both satisfied, under certain conditions. We note that the leading order gauge 
invariance of the truncated effective action has already been obtained by a different method in . 

The paper is organized as follows. In section (0 we discuss gauge dependence in the context of the 1PI 
effective action. In section ()3.1[1 we review the 2PI formalism. In section (|3.2|) we derive the 2PI Nielson 
identities. In section J3J we discuss the circumstances under which truncated 2PI equations can be expected 
to lead to gauge dependent results for physical quantities. In section J5J we look at the background field 
gauge and in section @ we present some conclusions. The connection between our result and that of is 
discussed in the appendix. 

Throughout this paper we use the compact notation of DeWitt A single latin index of the form 
{i, j, ■ ■ •} indicates the discrete group and Lorentz indices, and the continuous space-time variable. For 
example, a gauge field which would normally be written A®(x) becomes fa. Greek indices of the form 
{a, j3, • • •} indicate discrete group indices and space-time variables. For example, the Lorentz gauge 
condition, which would normally be written F a (x) — d^A a ^{x) becomes F a . The summation convention is 
used throughout, and is extended to include integration over continuous variables. 



2 The 1PI Effective Action 

2.1 Generalities 

We start from the partition function: 



Z 1PI [J] = Jvip Det \M aP \ exp [» (i + TtpA ] , I = S + S gf , (1) 

where S is the matter and gauge field action. The gauge fixing is set by the action S g f — ^F a F a . The 
matrix M a p is the Faddeev-Popov operator 



and the functions D%(ip) represent a complete and independent set of generators of the local gauge trans- 
formation. The only restriction on the gauge fixing condition is that the Faddeev-Popov matrix be 
invertible: we define the ghost propagator by 

MafiQ^ = —8 ai , (3) 

where the <5 Q7 in this equation is a Kroniker delta. 

It is usually more convenient to work with the effective action, instead of the partition function. We 
start from the generator for connected green functions: 

W[J] = -ilnZ . (4) 

The expectation value of the field is obtained from 

SW 

V = (<?) = jj-- (5) 

The effective action is defined by the Legendre transform 

T[fa = W[J] - J t fa ; (6) 
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Using J3J) we obtain the equations of motion, 

sr\ 



(7) 



2.2 Nielsen identities and gauge dependence 



The gauge dependence of the effective action can be explicitly calculated. We consider an infinitesimal 
change of the gauge condition: 

F a —> F a + 6F a (8) 

The goal is to calculate the change produced in the generating functional: 

SW = W[F a + 5F a ] - W[F a ] . (9) 

and to obtain the corresponding change in the effective action by Legendre transforming. The calculation 
can be done in a straightforward way by observing that the change in the action / produced by (JSJ can be 
canceled by a transformation 

f — ► ip + Sip , (iO) 
which amounts to a shift of integration variables. We define Sip through the equation, 

6(F[<p]) = F{ip + Sip) + SF(p + Sip) - F(ip) := . (11) 
Expanding to first order and introducing the notation F"{ip) = we have, 

F^Sipi = -SF a (ip) . (12) 

The unique solution of l|12|) is, 

Sip* = Di[ip]g%[ip]5FP[ip] . (13) 

Since this expression has the form of a gauge transformation, the gauge fixed action S + ^(F[ip\) 2 and the 
measure T>ip Det|M| will be invariant under © and p()ll3|l . Thus, the only contribution to 5T will come 
from the source term. We obtain, 

SW \ J=const. = W l F<X + 5F ^ - = MSP') 

8T\ T =SW\ 7 - Ji5(j) 1 

I J —const. I J —const. T 

54>i = (Sipi) + iJ k (Sip k ipi) ~ iJ k (Sip k )(/)i . (14) 
Combining we obtain the Nielsen identity for the 1PI effective action: 

^\ J=cons i=-X^,W-VW) (15) 

Thus we find that on the mass shell (defined by J; = — = 0) the effective action is gauge invariant. 

We note that the variation of the effective action in l|15[) is different from the expression found in Eq. 
(2.14) of Although both results describe the variation of the effective action caused by a change in the 
gauge condition, Eq. I|15|l is obtained by holding the source J constant, and Eq. (2.14) of |3j is obtained 
by holding the mean field cf> — </>[</, F] constant by varying J. It is easy to obtain the relationship between 
these quantities. We start with the generating function which we consider as a functional of the gauge fixing 
function and the source: W[F, J]. The variation of the generating function is given by, 



■SF+ — 

J— const. OJ 



■SJ 

F— const. 



S W\ J=const +4>SJ (16) 
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Subtracting <5 (</</>) from both sides and using the definition © we have, 

^ = ^L=con St .-^ ( 17 ) 

Considering the effective action as a function of the gauge fixing function and the expectation value of the 
field we have 



ST 

SF 



<p— const. 0(p 



F— const. 



= 5T\, . -J8<j) (18) 

\<p— const. T v ' 

Comparing fTTjl and (fl%)) we obtain, 

5T\, = 8W\ J = -TiiStp*) (19) 

\{p— const. \J=const. ' * ' ' x 1 

which can be compared with (|15f) . Throughout this paper we will consider variations obtained by holding 
sources constant. For the 2PI case in particular, this seems to be the more physical choice. The subscripts 
J = const will be dropped throughout. 

3 The 2PI Effective Action 
3.1 Generalities 

The 2PI action functional ^Oj can be defined using path integral methods following the technique used to 
construct the standard 1PI effective action. A bilocal source K is needed in addition to the standard local 
source J to define the partition function Z[J,K]. The generating function of the connected green function 
W[J, K] = —i\nZ[J,K] is also a two variable function: 

Z[J, K] = e iW ^ = J v^eiiSM+Ji^+^Wi) (20) 

The mean field <f>i and the connected two point function Gy are obtained from W[J, K] as 

6W S 2 W 

* = ^ = wr GV = ( ^ V> - ^ = • (21) 

Differentiating W[J, K] with respect to the bilocal source gives the following relation between the mean field 
and the connected two point function: 

SW 

+ 0V = 2f£- . (22) 
dKij 

The 2PI effective action functional is the Legendre transformation of W[J, K] with respect to J and K: 



T[<t>, G] = W[J, K] - jrf - -Katfp + G ij ) . (23) 



Using H21JI and (|22(l we obtain the following relations: 



W.G] t ±1. S T[cb,G} _1 

5G*i ~ 2 



Ji-Kirf; -^-±=--K l0 . (24) 



The effective action is usually written |1U| in a more convenient form which has a simple diagrammatical 
interpretation in terms of 2PI diagrams 

m G] = S [4>} + *i Tr {log (G- 1 ) + G (G 1 - G" 1 ) } + $[0, G]. (25) 
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where So is the free part of the action and Go is the bare two-point function (— iS 2 So[(p}/ StpSip^) 1 . The func- 
tional $[(/>, G] is the sum of all two-particle-irreducible (2PI) skeleton diagrams with bare vertices and dressed 
propagators. In the absence of sources, the difference between the inverse dressed and bare propagators is 
proportional to the one-particle irreducible self energy 

cr* -a? (26) 

This is the usual Schwinger-Dyson equation for the propagator. 

The above procedure can be generalized .11 to construct NPI effective actions. For the NPI effective 
action the skeleton diagrams are calculated using dressed propagators as well as dressed N-point proper 
vertices; the (N+l)-point vertex is bare. 

3.2 Off-Shell Gauge Dependence 

To study the gauge dependence of the 2PI effective action we derive the corresponding Nielsen identities. 
For simplicity, we consider a pure Yang- Mills theory. The partition function with a gauge fixing term F has 
the form, 



Z[J,K] 



_ e iW[J,K] 



T>(p Det 



5F a \ip] . 



Stp t 



< xp / ( Sym + \{F[ V ]f + J*<pi + ^;,A 'V,) . (27) 



We calculate the variation of the generating functional W[J, K] under the transformations (JSJ and I|1UI13[I . 
As in the case of the 1PI effective action, the measure and the gauge fixed action are invariant which means 
that the only non-zero contribution comes from the source terms. The sources themselves (J, K) are viewed 
as the independent variables, and are held constant. We obtain, 

SW[J, K] = W[F + SF] - W[F] = Ji(V) + ~Y [(<fyV> + (V'V)] (28) 
The variation of the effective action is given by 

ST = SW - JiS4> - ^6(<p m ) . (29) 
Using the functional definition of expectation values, we obtain 

SW = (T) 

Sfc = (6<pi) + i(ifiiT) - i4>i(T) 

5{ipnpj) = (ifiiSipj + Scpupj) + i(ipiipjT) - i(tpnpj)(T) (30) 

where we have defined, 

T = JiSip* + ^Kijtf&tpi + <W) (31) 
Combining these results and using l|24|l we find: 

T--'((£* + 3^W'« , + Vrt) (^(' + §j8>')) , (32) 

where = ipi — and Gy = — Gij . Using (|24|l we see that the full 2PI effective action is gauge invariant 
on-shell. Of course, this conclusion is obvious, since the 2PI effective action calculated to all orders is exactly 
equivalent to the IPI effective action. One way to see this point is to note that the bilocal source K does 
not play an active role in (|2~TJl . We could set K to zero before differentiating, which shows explicitly that 
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the 2PI Gij is the same as the 1PI G^. 



Note that the variation of the 2PI effective action at constant mean field and propagator can be obtained 
following the method that was used for the 1PI effective action. We find: 

( ST ) G}=const. = ( SW ) I {J, K}— const. ' 

(33) 



4 Truncation 

All of the preceding calculations are valid when one works with the full effective action (1PI or 2PI). In 
practice, of course, we never calculate the full effective action: we use an expansion scheme and truncate 
at some finite order. When we work with a truncated 2P1 effective action, the non-perturbative nature of 
the 2PI formalism gives rise to problems with gauge invariance. In fact, it is easy to see in advance that 
this problem will occur. The 2PI formalism involves the use of resummed propagators. When calculating 
at finite orders, one effectively resums a specific class of topologies. Since the Ward identities involve the 
cancellation of contributions from different topologies, we expect that a resummation that involves only 
one type of topology will give rise to violations of the Ward identities. In this paper we study the gauge 
invariance of the effective action. The issue of the Ward Idenities will be left for a future publication |12j . 

We show below that both the IPI and 2PI effective actions are gauge invariant on shell to arbitrary 
order in any self-consistent expansion scheme. 



4.1 IPI 

First, consider calculating the full IPI effective action, i.e. to all orders in the expansion parameter. The 
mass shell condition is obtained from (JTJ with the source set to zero: T.i = 0. Substituting into l|15|l we have 
STipj = which tells us that the full IPI effective action is gauge invariant on-shell. We discuss below the 
truncation of the full IPI effective action. For definiteness, we consider a loop expansion. Calculating up to 
L loops gives the truncated effective action. The full effective action is the sum of the truncated effective 
action and the remainder: 



r = r L + r ex 



(34) 



where is calculated up to g 2L 2 and T ex ~ g 2L . When using a truncated effective action, the on-shell 
condition is replaced by an approximate on-shell condition determined from the truncated effective action 
by, 



= 



(35) 



Using this approximate on-shell condition we have 



8T 



r ex ~ g 



2L 



(36) 



Substituting into (|f 5fl we obtain, 



S(T L + T ex ) ~ g 



AL 



(37) 



which gives 



5T L ~ ST ex + 0(g 4L ) ~ T ex + 0(g 4L ) ~ g 2L 



(38) 



Thus we obtain, 



6T L ~ g 2 T L 



(39) 



G 



which shows that the gauge dependence of the on-shell effective action always occurs at higher order than 
the order of truncation 1 

It is important to note that the above formal analysis implicitly assumes that the solution of the 
truncated equation of motion Eq. I|35JI remains within the range of validity of the approximation implied by 
the truncation. In other words, it is possible that T ex , evaluated at the aproximate solution cf>o> contains 
terms of the same order as T^. In this case, the truncated effective action, and all physical quantities 
derived from it, can be gauge dependent. Such gauge dependence has been observed many times in the 
literature (the Coleman- Weinberg mechanism |6I7| . self-consistent dimensional reduction 8 , and the one 
loop plasmon damping rate |2I.''>| ). The gauge dependence in all of these cases is a manifestation of an 
inconsistent approximation scheme, since the gauge dependence identities guarantee that physical quantities 
will be gauge independent when calculated within a self-consistent perturbative procedure. 2 



4.2 2PI 

Now consider the full 2PI effective action (calculated to all orders). In this case, the conclusion we draw 
from 1132[| is the same as the conclusion from <|I5|) : on the mass shell, defined by ST/Scf) 1 — ST/SG^ — 0, the 
full 2PI effective action is gauge invariant. As mentioned earlier, this conclusion is obvious, since the 2PI 
effective action calculated to all orders is exactly equivalent to the 1PI effective action. Problems arise when 
we try to work with a truncated effective action. As in the 1PI case, we write the full effective action as the 
sum of the truncated effective action and the remainder as in 134|) . The approximate on-shell condition is 
determined from the truncated effective action by, 

Now let's consider performing a perturbative expansion of the type described in section l|4.I|) on i|32|l . 
Using H4l)f) we have 



ST , 6T e 



,2L 



(41) 



Note that in (|4I \ we have assumed that the effective action T ex and its derivative are of the same order in 
the expansion. In fact, it is not clear that this is the case, since the variable that we are differentiating with 
respect to is a non-perturbative quantity, however, it seems reasonable that differentiation will not increase 
the order of T ex . We also note that, as in the 1PI case, it is essential that the perturbative procedure is 
self-consistent, or that the quantity that has been dropped (5T ex {({P L , G°)) is higher order in g than the 
quantity that is kept (STl^'I, Gf)). Because of the inherent non-perturbative nature of the 2PI formalism, 
this possible mixing of orders is more difficult to monitor and control than in the usual IPI case. If we assume 
for the moment that the integrity of the loop expansion is not fundamentally violated, and substitute into 
(J2U, we obtain as before, 

6{T L +r ex )~ g 4L (42) 

which gives 

ST L ~ ST ex + 0(g iL ) ~ T ex + 0(g iL ) ~ g 2L (43) 

lr The variation of the effective action at constant can be evaluated using Eq.OHJ- In this case the variation of the effective 
action is linear in the source. Using the approximate on shell condition we obtain 

i^) <b=const. ~ TT ,n ~ 3 

d<p 1 4>l 

which shows that the gauge variation of total effective action is more weakly suppressed when the expectation value of the field 
is held constant. Note however that the result for the truncated effective action 1391 still holds. 

2 Note that a self-consistent perturbative procedure does not necessarily guarantee accuracy. A self-consistent perturbative 
expansion is one in which an expansion can be carried out in some parameter (such as h or a coupling constant g), in such a 
way that there is no mixing of orders in the expansion. If this self-consistent perturbative expansion exists, the existence of 
the gauge dependence identities ensures gauge independence of the on-shell effective action order by order, irrespective of the 
relative magnitude of each term in the expansion. 
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Thus we obtain the same result as in the 1PI case: 

ST L ~ g 2 T L (44) 

which shows that the gauge dependence of the on-shcll effective action always formally occurs at higher order 
than the order of truncation, just as for the 1PI case. 



5 The background gauge 



In this section we review the background field gauge technique. A more detailed description can be found in 
|13j. For simplicity, we consider a pure Yang-Mills theory. We write the gauge field as the sum of two pieces: 



gauge field = + Q p 



(45) 



where is a background gauge field and is a quantum gauge field which will be the variable of integration 
in the functional integral. We add a gauge fixing term (the 'background field gauge') that breaks gauge 
invariance in Q M but not A^. Also, we couple only Q a to sources. This procedure allows us to calculate 
quantum corrrections and keep explicit gauge invariance in the background field variable. 
We define our notation as follows: 

[T a ,T b ]=if abc T c ; if abc = (T b ) ac 



Dl c = 5 ac d^ 



D 



gfabcA b 



M - d n - igA°T a =8^- igA^ 
-igF av = [D^Du] 

= df+A" — d v A a ^ + gfabeApAl 

Where T a are the generators of the Lie algebra of the gauge group. 
The partition function is defined as, 



(46) 



Z[J,K,A] 



[dQ]det 



5G a 



Sa b 



exp 



\J d 4 x (c YM (A + Q) 



2£ 



(G a Y + JZ(x)QZ(x] 



} , d 4 x I d 4 yK^(x,y)Q^(x)QUy) 



(47) 



We define the generating functional for connected diagrams, and the effective action, in the usual way: 
W[J,K,A] = -ilnZ[J,K,A] 



T[Q, G, A] = W[J, K, A] 



d A xd A y^K^Q» +&") 



sw 



Guv 



S 2 W 



8{iJ u )5{iJ v ) 



We use the background field gauge condition: 



G a = M Q£ + f abc A b Q£ = (D^) a 



(48) 



(49) 



Note that the propagator G^ u should not be confused with the gauge fixing functional G a . We consider the 
transformation: 



AO, 



A a 



5A«: 



D u a b 



SQu = fa 



be I 



(50) 



Note that the transformation on A has the form of a gauge transformation. The transformation on Q is just 
a shift of the integration variable. It is straightforward to see that Q5O0 gives, 



6{A«+Ql) = Dfa c] D a = 8 u 
S(G a ) = f abc G b a c 



ig(A^ + Q„) 



(51) 
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which means that both the Yang-Mills Lagrangian and the gauge fixing term remain invariant under l|5U|) . 

Combining these results we see that when the transformation (|50() is performed on the generating 
functional or the effective action, the change that is produced comes only from the source terms. Note that 
this situation is exactly analogous to the situation we had previously. In sections J5J and © we shifted the 
gauge fixing function and performed a simultaneous shift of the integration variable so that the changes to 
the generating function and the effective action came only from the source terms. In this case, the shift of 
the gauge fixing function is generated by a shift in the background field, as shown in l|51() . The calculation 
of the Nielsen identity follows the same procedure as before. 

We can calculate the Nielsen identity for T[Q,G,A] by following the procedure in section (|3.2|) . We 
perform the change of variables (|5t)(l and obtain, 

ST = T[Q, G,A + 5 A] - r[Q, G, A] 

= -•((^«' + i&wa +a « i aO(^ + lf^)) (52) 

where £,\= Q l — Q 1 [A], and G l] A — Ca^a ~ Gi 3 [j4]. As before, we find that the full effective action is gauge 
invariant on shell, and the truncated effective action is gauge invariant to leading order. 

6 Conclusions 

We have derived the gauge fixing identities for the 2PI effective action, valid for any gauge fixing function 
with an invertible Faddeev-Popov matrix. These identities were first derived by Arrizabalaga and Smit [Sj 
using a different formalism. As expected, these identities prove that the 2PI effective action is invariant 
under infinitesimal gauge variations, on shell, to arbitrary order in any self-consistent expansion scheme. 
We have also considered the background field gauge and shown that the effective action is invariant under 
infinitesimal shifts on the background field, on shell, to arbitrary order in any expansion scheme. 

We note that the derivations of the gauge fixing identities for the 1PI and 2PI formalisms are virtually 
identical from a mathematical point of view. This simularity is unexpected, in light of the fact that we 
expect gauge invarience problems of a completely different nature to arise in the 2PI theory. The 2PI theory 
is inherently non-perturbative, and involves not just the mean field 4>i but also the two point function Gij, 
which are a priori independent, and must be solved for simultaneously. As a consequence, the 2PI theory 
preferentially resums specific topologies, a procedure that we expect will lead to violations of Ward identities. 
Of course, since Ward identities reflect the quantum symmetries of the Green functions, such violations are 
due to truncations and would be absent if exact calculations were possible. 

One possible problem with the truncated 2PI effective action is the fact that the perturbative procedure 
is not necessarily self-consistent, in the sense discussed in sections 4.1 and 4.2. It is also possible that there is 
a fundamental problem associated with taking the Legendre transform of a truncated theory. In general, the 
Legendre transform is defined functionally as the transformation that converts the untruncated generating 
functional to the untruncated effective action. The 1PI theory is inherently perturbative and thus there 
is no problem with the Legendre transform for truncated versions of the theory. However, in the case of 
the 2PI theory, it is possible that problems arise when transforming the truncated theory. In this case, 
the minimum of the effective potential would not necessarily correspond to the expectation value of the 
generating functional, and the interpretation of the Ward identities obtained from the effective action would 
not be straightforward. 



Appendix A: BRST transformation 

Our result agrees with the result of [S] which was obtained using the BRST method. The basic strategy of 
the BRST method is to exploit the fact that the gauge-fixed theory still possesses a global symmetry called 
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the BRST symmetry. This symmetry is made explicit by using a representation of the partition function 
that is different from J57J. The gauge fixing term has the form, 

Sgf = J d A x (^-c a M al3 C(3 + B a V a - ]pcB a B^j (53) 

where c a and c a are the ghost fields, B a is the auxiliary field, and V a is the gauge fixing condition. Integration 
over the ghost fields produces the determinant of the Faddeev-Popov matrix, as in (|27[) . and integration over 
the auxiliary field produces a gauge fixing term of the form ^V a V a . Comparing with l|27|) we have, 

V a = ^F a (54) 

To make further progress one must specialize to the covariant gauge: V a — ► d^A'tix). The gauge fixed action 
(Sym + Sgf) is invariant under the BRST transformation: SbrstA® = e(D M c) a , Sbrstc 11 = ieg(T a c a )(TbC b ), 
Sbrstc 11 = —eB a , and SbrstB 11 = 0, where e is an infinitesimal global anti-commuting parameter. One 
obtains; 



Sgf = Qbrst J <Fx [^Xc a B a - c a V a j := Qbrst* (55) 

where Qbrst is the BRST nilpotent charge defined as Sbrst^P — zQbrst^P- Calculating the variation of 
the effective action under the BRST transformation one obtains, 

STbbst - \ Uqbrst + §jG-) 2 ) , (56) 



with 



= - / d 4 x ( c a SV a [A] - X -b X c a B a J ( r,7 ) 



It is straightforward to see that this result is equivalent to ours. We take V a — ^fxFa which gives 

5F=^-5V-\-^-V (58) 
We integrate over the Gaussian B field, and over the ghost fields using 

J VcDcciCje-^ Mc) = [M-%detM . (59) 
We find that the variation STbrst is proportional to our effective action variation <5I\ 
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